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I. INTRODUCTION 


Black hole hair has been a topic of discussion in the literature for over forty years now ever 
since Israel and Carter’s original uniqueness (“No-hair”) theorems for static, asymptotically 
flat black holesi”- which classihed all black holes by their mass m and electric charge e. 
This became a subject of active research in the 1980’s due to the discovery by Bartnik and 
McKinnon^ and Bizon^ that these theorems could be violated by considering Einstein-Yang- 
Mills (EYM) theory, which imposes extra gauge symmetries on space-time - these manifest 
as gauge helds coupled to the gravitational held. The spirit of the conjecture was preserved 
in a way, since the known solutions were still specihed by a relatively small number of 
parameters, and for asymptotically hat space the solutions were found to be sparse, existing 
with only discrete values of boundary conditions'^-. In addition these are unstable, since 
it is known that the number of unstable modes of these solutions is proportional to the 
number of zeroes of the associated gauge held function (“nodes”), with 2n unstable modes 
for gauge functions with n nodes, and in addition, that solutions with A = 0 have at least 
one node^^— . Much work exists in the literature for hairy black holes in the case A = 0 
involving a wide variety of special cases, including non-spherically symmetric geometries and 
higher-dimensional analogues^^^^— . 

The dihering geometry bestowed by a negative cosmological constant yields interesting 
results. For one, we hnd that the solution space is a lot more rich and abundant, with 
solutions existing in open sets of the initial conditions as opposed to discrete families^^— ; 
and we may hnd ‘nodeless’ solutions (i.e. which possess gauge held functions with no zeroes), 
which we have said are of importance to stability. When we take A < 0, we hnd solutions 
where the gauge held function uj has no zeroes and which are stable under linear—^— and 
non-spherically symmetrio^*^ perturbations, provided we take the limit |A| oo. 

The case of hairy particle-like or ‘soliton’ solutions has also been considered. The non¬ 
existence of gravitational^^ or pure Yang-Mills^ solitons in asymptotically hat space preceded 
the surprise discovery of EYM solitons in su(2) asymptotically anti-de Sitter (adS) spaced. 
We note that these are only stable for A < OiS, so naturally this inspired generalisations 
of their work in asymptotically adS spaced. Indeed, solitons have been a subject of much 
interest in a range of special cases involving existence and stability^^— . 
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These ideas have been extended in many ways. One obvious idea is to extend the 
gauge group from su(2) into su{N), which is something that that has been worked on 
previouslyi^'^^TiS^. Another extension is to move away from the assumption of spherical 
symmetry to considering solutions on manifolds of alternative topology^^^— . Notable here is 
the work of van der Bij and Radi>^, who considered manifolds with isometry groups based on 
foliations of space-time by surfaces of constant Gaussian curvature: these are parameterised 
by an integer k = 1,0, —1, the sign of the curvature; and the foliated surfaces of the space 
have a spherical {k = 1), planar [k = 0) or hyperbolic {k = —1) isometry group. That work 
inspired a previous work by the suthor— , and these ideas will appear here. We note early 
on that while black hole solutions are possible for all 3 values of k, solitons are only possible 
for fc = 1, the spherical case: this is because for k = 0,-1, the Riemann curvature scalar R 
blows up at the origin, and therefore there is no such thing as a globally regular solution for 
k 1 throughout the range [0, oo). 

Another subject of interest has been in so-called dyonic solutions, which have a non-zero 
electric sector of the gauge potential, unlike previous work which has primarily concerned 
‘purely magnetic’ solutions. Non-existence of genuinely non-Abelian regular monopoles and 
dyons (the dyonic analogue of solitons)^ was proven for flat space su(2) EYM— , but then 
later black hole solutions and stable monopole and dyon solutions were found in asymp¬ 
totically adS 5u(2) EYM— In addition, dyons and dyonic black holes have been found 
in spaces with axial symmetry^>^. Notably, Nolan and Winstanley^ recently proved the 
existence of dyons and dyonic black holes in four-dimensional su(2) EYM theory with A < 0. 
It is our intention to extend this work in four-dimensional dyonic solutions to consider an 
5u{N) gauge field, and for topological black holes of the kind considered by van der Bij and 
Radii^. 

In this paper, we prove the existence of four-dimensional, topological, dyonic black hole 
and soliton solutions to su{N) EYM theory with A < 0. In Section [TTl we review the 
features of the model in question, and derive the field equations in this case. We also prove 
the existence in Section III Bl of several trivial solutions to the equations, whose existence 
has been proven elsewhere. In Section HTTl we prove several propositions concerning the local 
existence and analyticity of regular solutions near the spatial boundaries of the solutions, i.e. 
near r = 0, r = Vh and as r —)■ oo. Then in Section HVl we prove some propositions concerning 
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the global behaviour of solutions: notably, that solutions which start regularly near the event 
horizon for black holes (or the origin in the case of solitons) can be integrated arbitrarily far 
out into the asymptotic regime, as long as the metric function /i > 0; and that solutions in 
the asymptotic regime will remain regular as r —)■ oo. We use these propositions in Section 
El to show we may ‘stitch’ locally existing solutions together into global solutions - this is 
a fairly standard ‘shooting’ argument that has a well-established usage in the literature^i^. 
Essentially, we prove that solutions to the dyonic held equations exist in open sets. This 
allows us to use some trivial solutions and some previously established existing solutions to 
deduce the existence of non-trivial solutions to the held equations in various regimes. We 
also prove in Section IV Bl that we can hnd solutions in the limit |A| —)• cxo for arbitrary initial 
values for one gauge held and initial values in a neighbourhood of zero for the other. We 
hnish in Section |VT] by presenting our main results and hnal conclusions. 

II. 4D TOPOLOGICAL su(iV) EINSTEIN-YANG-MILLS THEORY 

In this section we give the action and ansatze for topological su{N) EYM theory for 
A < 0, give a brief derivation of the dyonic held equations, and list some trivial solutions in 
this case. 

A. Ansatze and field eqnations 

The action we shall use for four-dimensional su{N) EYM theory with a negative cosmo¬ 
logical constant is 

Seym = - j d^x^-detg[R - 2A - TrE^^FH- (1) 

Here, R is the Ricci scalar of the geometry; we take A, the cosmological constant, to be less 
than zero; and the anti-symmetric held strength tensor is given by 

Ffiu = d^Ay — dyAfj_ + q [A^, Ay]. (2) 

In all that follows we use units in which AttG = c = q = 1, and the 4D metric we use has 
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signature (—, +, +, +). Varying the action ([T]) gives the held equations 



( 3 ) 


where the Yang-Mills stress-energy tensor is 



(4) 


Note that the square bracket in ([3]) is the Lie algebra commutator, and ‘Tr’ in ([H 0]) is the 


Lie algebra trace. Throughout this work we have employed the usual Einstein summation 
convention, unless otherwise stated. 

In this paper we focus on static, topological black hole and soliton solutions of the held 
equations ([3]), specihcally for spaces regularly foliated by 2D (spacelike) hypersurfaces of 
constant Gaussian curvature which can be indexed by /c G { — 1,0,1} for positive, zero, 
and negatively curved spaces^. Hence, we may write the metric in standard Schwarzschild 
co-ordinates as 


ds^ = + /i ^dr^ + r‘^dVt\ 


(5) 


where /i and S depend on r alone. We may write the function /i as 



( 6 ) 


noting that the constant ^ is the adS radius of curvature dehned by 



possible because we are only interested in A < 0. We will refer to both forms of the 
cosmological constant in this work. The angular part of the metric, dVL\, is given by 


dnl = de^ + fl{e)dct?- 


( 8 ) 
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where we recall that if we let 


fk{0) = 


sin^ for k 
6 for k 
sinh0 for k 


1 

0 

-1 


(9) 


then our metric ([5]) is endowed with the topology appropriate to each case. 

In previous work— we only considered ‘purely magnetic’ solutions for topological su{N) 
black holes. Here, we wish to include the electric part as well, meaning that we use the full 
gauge potential, which in our case is 


A = A^dx'^ = Adt + I3dr + ^(C- C^)d0 - ^ 


(C + C^)M0) + D 


dfkie) 

dO 


d(j), ( 10 ) 


where A, B, C and D are all {N x N) matrices, is the Hermitian conjugate of C, and 
fk{0) is given above as ([2]). The matrices A and B are given by 


[^]jj = [^]jj = 


( 11 ) 


for 2N functions aj{r) and Bj{r) {j = 1, constrained by 

N N 

5^a,(r)=5^B,(r)=0. (12) 

i=i j=i 

Hence these represent 2N — 2 independent functions. The matrix C (which also depends 
solely on r) is upper-triangular, with non-zero entries only immediately above the diagonal, 

i.e.: 

|CL+i=(|Cy+,,,)* = c,(r)e«>'’') (13) 

(where * is the complex conjugate) for 2N — 2 functions u:j{r) and gj{r) {j = 1,..., iV — 1). 
The matrix D is constant and diagonal: 


[B]jj — N + 1 - 2j. 


(14) 


We hnd that one of the Yang-Mills equations gives us Bj + gj = 0, hence we let gj = —Bj. 
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Finally we may exploit a remaining gauge freedom to set Bj = gj = Therefore we have 
2N — 2 independent gauge held functions altogether: aj(r), which we shall call the Electric 
Gauge Functions (EGFs), and u}j{r), which we shall call the Magnetic Gauge Functions 
(MGFs). For further details on this choice of potential, note that the precise form has been 
derived in detail^^ following a method due to Kiinzle^^. 


Substituting the metric ([5]) and the gauge potential flTOD into the held equations 
using the Bianchi identities) gives us the two Einstein equations 

^2 A 

S' _ C 2G 
'S ~ 2iFS^r ^ V’ 

where for convenience we have dehned new quantities 

N 

i=i 

N-l 


7V-1 


i=i 


V 




(and 

(15) 

(16) 


(17) 

(18) 

(19) 

( 20 ) 


i=i 


and 21S — 2 independent Yang-Mills equations 




(S' 2 \ 1 

[j --J ^ - «i+i) - («i-i - «i)) ’ 




fi' 

b p 






-Uj {oij 


WiU 




lir^ 


( 21 ) 

( 22 ) 


where we have dehned 


Wj = k- a;| -1- - 


(23) 
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Note that a prime ' stands for d/dr, and we define ao = aN+i = ojq = = 0. Also, 

it may be noted that there are N eqnations in aj, bnt only — 1 of the eqnations are 
independent dne to the degree of freedom given to ns by (IT^ . If we take the above held 
eqnations (na [ig ED 122]) and let = 2 and fc = 1, we can verify that they rednce to the 
correct limit for dyonic spherically symmetric 5u(2) EYM held eqnations^; and if we let 
aj{r) = 0 we recover pnrely magnetic topological su{N) EYM eqnations^^. 


As may be expected from previons cases^^'^'^, the held eqnations flT^ [TBl 1211 122]) are 
invariant nnder varions transformations. They are invariant nnder the following two trans¬ 
formations separately: 

^ (Vj), j ^ N -j (Vj); (24) 

they are invariant nnder 

OJj I—^ 

for each Uj separately; and hnally they are also invariant nnder the simnltaneons transfor¬ 
mations 

S{r) !-)■ AS'(r), (^j{r) i-)- Xaj{r), (26) 

which correspond to applying a time-rescaling 1 1 —)■ X~^t to the metric (|5|) and the potential 

oa. 


Finally, we re-express the eqnations in another form which will prove nsefnl later. We 
make the variable change 

~ ~ ®j+i- (27) 

This leaves the forms of flT^ [TB]) nnchanged, thongh it does alter ( flTS]) : 

N-l 

C s (28) 


It also alters rj bnt we deal with that as it comes np as the transformation is qnite complicated. 



The Yang-Mills equations are altered to become 


S'; 




SS' 2\ 1 

(S' 1 ..,2 

is fi Y 4fi2S2 ^ J fir2 


(29) 

(30) 


with 

= 2ujSj - - ^‘^^^^Sj+i. (31) 

We should mention that since fl271) is an elementary linear transform, it turns out to have 
an elementary inverse: 

-1 ^"V h\ 

^ + Sk. (32) 

fc=l k=j ^ ' 

It can be noticed that this equation has been stated before^^ as a way of expressing the 
functions ctj, and this connection is to do with the fact that, in vector form, we can express 
the transform (l27)l as £ = Ta, where T is the {N — 1) x N matrix 


T = 


/l -1 0 0 ... 0 0 \ 
0 1 -1 0 ... 0 0 
0 0 1 -1 ... 0 0 


\0 0 0 0 ... 1 - 1 / 


(33) 


whose rows represent a choice of simple roots for su(iV)^, and expressing the equations thus 
will reveal the structure of the Cartan matrix of 5u(Y) in the electric gauge equations - as 
we shall see, the form of (l3T|l simplihes the analysis considerably. 


We lastly note that we can also use the Einstein equations to express fl2^ in the following 
useful form: 


= -2(m-rP + -- — 






(34) 
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B. Trivial solutions 


We may find the following list of trivial solntions, the existence of most of which has been 
investigated in previons work. 


1. If we set aj{r) = 0, we recover the pnrely magnetic topological solntions which we 
previonsly investigated^ (thns they exist for all k). We note that in this case, eqnation 
(fT^ deconples from the others. 

2. If we set aj{r) = 0 and ujj{r) = \/j{N — j), we obtain the Schwarzschild-anti-de Sitter 
(SadS) solntion where m{r) = M and S{r) is a constant which we nsnally scale to 1. 
(This solntion only exists for k = 1.) If M = 0, we recover pnre adS space as a solntion. 


3. If we set aj(r) = ojj{r) = 0, we obtain the magnetically charged Reissner-Nordstrom 
anti-de Sitter (RNadS) black hole solntion. Here, S'(r) is again a constant bnt 


m{r) = M 


Qm 


(35) 


and the magnetic charge Qm is given by 

k^N{N + 1){N -1) 
Qm — - - - 


(36) 


This solntion exists for all k. 
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If we instead set u}j{r) = 0 and aj{r) = — (with constants aj), we get (for all k) an 
su(A^) analogy to the su(2) Abelian RNadS black hole^ with S{r) = 1, 


m{r) = M 


Qe + Qm _ 
^ ’ 


( 37 ) 


where magnetic charge Qm is given in fl36ll . and total electric charge Qe is 






in which the individual charges aj are constrained by 


N 

= ( 39 ) 

i=i 

We note that apart from fl3^ there appear to be no a priori constraints on the indi¬ 
vidual OjS, though one possible choice is Oj = Q{N -|- 1 — 2j) for some constant Q: this 
represents the 5u(2) embedding of the solution (see (HIT) ). Interestingly, this choice 
would give us 

0 . ^ Q^NiN 

0 

so that Q^Qm = k'^QE- The reason this is of interest is that for fc = 0 we can evidently 
hnd solutions with a zero magnetic charge and non-zero electric charge, which we know 
is impossible in the spherical case^. Due to the fact that ooj = 0, we might expect that 
this solution would not be stable under linear perturbations (in analogy to the su(2) 
results^^— ; though very recent results we have obtained suggest that for k ^ 1 this 
may not be the case after all^. Finally, we note that this solution has not yet been 
investigated, though it is worth noting for future reference. 


Finally, we can also obtain dyonic su(2) embedded solutions with the following rescaling: 


Proposition 1 Any solution which satisfies the3u{2) dyonic field equation^ can he rescaled 
and embedded as an su{N) dyonic EYM solution (which satisfies KTR \TU. [HI 


Proof We begin with the held equations flTSl[T6lIHll22|) . We rescale them with the following 
dehnitions: 


A^^iiV(iV-l)(iV + l), 
^ \/j{N - j)u: Vj, 
aj^{N + 1- 2j)a Vj, 


r = XjsfR, i = 
m = XnAi, 

S = XnS. 
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(41) 





This rescaling leads to the following equations: 


dm /da\^ u'^a^ /da;\^ —/c)^ 

^ 2i?2 

1 dS 2u‘^a‘^ 2 /da;\^ 

^d^ “ /i2^2/2 + R ’ 

d^a /1 d^ 2 \ da 2a;^a 
d^ ^ ~r)^^ fiR^ ’ 

d^ca /1 d/i 1 d^\ do; a^o; a;(a;^ — 

iff “ ~ 1 fidS ’*' IdR I dS ~ ^252 ^Ir2 


( 42 ) 


It can be checked that if we let /c = 1, these are exactly the same equations as the dyonic 
spherically symmetric su(2) held equations for which existence of solutions has been proven?^; 
and if we let a = 0 then we recover the original topological su(2) equations^. We also note 
that the tracelessness condition fll2p constraining the aj is still satished.D 


In this paper we investigate the existence of solutions to the held equations (1T5| ITHl [2ll 
[22I) which do not appear in the above list - i.e. genuine non-trivial solutions. 


C. Boundary conditions 

As we are interested in both black hole and soliton solutions, the boundary points we must 
examine are the origin r = 0 for solitons, the event horizon r = Vh for black holes, and for 
both solutions, the limit r —>■ 00 . The held equations are singular at all of these points, so we 
shall now assume appropriate power series solutions regular near the boundary in question, 
substitute them into the held equations, and also use some requirements of physicality (e.g. 
expected asymptotic behaviour, regularity) to determine the local power series expansions - 
in particular, we are interested in how many independent parameters are required to specify 
any solution. 
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1. Origin r = 0 


At the origin, we simply use the co-ordinate r. The following power series forms are 
assumed: 


m(r) = mo -|- mir -|- m 2 r^ -|- O(r^), 
^(r) = ^o + ^ir + ^2r' + 0(r3), 
u:j{r) = uijfi + Uj^ir -f- -h 0{r^), 

£j{r) = 8 jfl + Sj^ir + 8 j^ 2 r‘^ + O(r^). 


We then substitute these expansions into the held equations flT^ [TBl 122112^ . To avoid a 
singularity in the metric and the held equations, it must be the case that Sq is non-zero. 


mo = mi = m 2 = Si = Sj^ = ujj^i = 0, and 

= \/i(N -j). 


(44) 


for all j G {1,..., — 1}. The metric function expansions are not difficult to compute. The 

Yang-Mills sector is vastly more complicated and involves solving a matrix equation. We 
shall look into this in much more detail later in Section UlIAj for now we shall merely present 
the results. We dehne vectors uj = {uji,uj 2 , ...^ojjsi-iY' and S = {£ 1 , 82 , ■ When we 

substitute the conditions into the held equations, the analysis implies that the expansion 
is best done in terms of the N — 1 eigenvectors of the matrix in question fl73|) . so that in 
order to obtain all the independent parameters, we must expand n; up to order and 8 up 
to order r^~^. The expansions nearby the origin fjTdll thus become (in components) 


m(r) = m^r^ + 0(r'^), 

^(r) = So + S2r^ + 0{r^), 


N-l 


=VjiN -j) 1 + 5^ f3kix)vy+^ j + 0(r^+^), 

\ k=l 

N-l 

YY) = X] h{x)vlr^ + O(r^). 


( 45 ) 


k=l 
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We will define the vectors vl properly later; just now all we need state is that /3k{x), 9k{x) are 
determined entirely by the 2N — 2 parameters ^^(0) = and 0^(0) = Ok] and the constants 
m 3 and S 2 are determined entirely by Sq, (or Sj^i) and Uj ^2 by the held equations as 
follows. 

/ iV 7V-1 \ 

""3 = - J)^ll 

° Vi=i i=i / 

4 1 ^ / 2 

+ Y + 6 J>i,2- \/(j-i)(^-i + iVi-i,2) , (46) 

j=i i=i 

.. N-l N-l 

° i=i i=i 

The value of Sq is hxed by the requirement that S' —)■ 1 as r —)■ 00 , therefore these solutions 
are specihed entirely by 2N — 2 parameters. 


2. Event horizon r = 


At the event horizon, the picture is much clearer. Firstly we have /i(r/i) = 0. This means 
we immediately have 


kvh 


rrih = 


H- 

2£2 


(47) 


It also means that (12111 will be singular at r = r/j, unless Ojirh) = 0 for all j. Substituting 
these two facts into the held equations flTHl lT6l l2ll l22ll produces the following regular Taylor 
expansions in the held variables at the event horizon: 


m(r) = mh + — Vh) + 0{r — ThY, 

S(r) = Sh + S'^{r - Th) + 0{r - Vhf, 
aj{r) = - Vh) + 0{r - Vhf, 

ujj{r) = Uj^h + uj'j f^ir - Th) + 0{r - Vhf] 
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where the expressions for m'f^, are given in terms of Uj^h and Sh- 



Therefore and are entirely determined, and we also have 

, ,, , k 2ml 3r/, , , 

k‘h = k- i'^h) =-h > 0, (50) 

Th Th r 

where the last inequality is due to the requirement of a regular event horizon. Thus, fixing 
r/j and £, and hxing Sh with the requirement that S' —)■ 1 as r —)■ cxo, the solutions near the 
event horizon can be entirely specihed by 2N — 2 parameters. 

We should also finally note a weak constraint on ujj and aj near the horizon, given by 
requiring a non-extremal event horizon: setting fih = C(j,h = 0 in the held equation for m' 
gives us 

+ 2F(rJ <k + !^, (51) 

as directly related to the tangential pressure^. This also allows us to dehne a minimum 
event horizon radius for k = —1, 

rl>j (2mi + 1), (52) 

where it can be seen that the right-hand side of the inequality is clearly positive. Finally, 
we note that it gives us a minimum bound on |A| for k = —1, in analogy to the su(2) case^: 

|A|>^(l + 2P(r.) + :h|i)). (53) 
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3. Infinity 


As r —)■ oo, we wish the solution to approach ‘topological adS’ - that is, we wish m ^ M 


2 

(some constant) and S' —?• 1. We also note that in this regime, y{r) —)■ hence using fITB]) 


shows that S'{r) ~ 0{r ®). Therefore, the expansions in the asymptotic region are: 


m(r) = M + mir ^ + 0 (r ^), 
S{r) = l + Sir-^ + 0{r-^), 

aj{r) = + 0(r“2), 

Uj{r) = + 0(r "2), 


(54) 


where mi and S 4 are given in terms of aj^oo, dj, oJj^oo and cf 



(55) 


Note that it is easy to substitute aj for £j here using fl^ or ([32]), if we need to obtain the 
expansions in terms either variable. No conditions are placed on the constants ^j,ooi 
Cj, dj or M, hence (hxing Vh and £) we get a 4iV — 3 parameter family of solutions. 


III. LOCAL EXISTENCE OF SOLUTIONS NEAR BOUNDARIES 

We now begin by proving that solutions to flTSl [T 6 l [211122]) exist in some neighbourhood 
of the boundary points r = 0, r = Vh and r —)■ 00 , and that these solutions are analytic in 
their boundary conditions in some sufficiently small open neighbourhood of the parameter 
space. To do this we use a well-established a theorem of differential equations, which has a 
rich history of use in the literature^'^'^. We begin by stating the theorem. 
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Theorem 2 — Consider a system of differential equations for n + m functions a = 
(ai, 02 ,, a„) and b = ( 6 i, 62 ,..., bm) of the form 

= xP^fi{x,a,h), 

x^ = -Xih + x'^^gi{x, a, b) 
dx 

with constants Aj > 0 and integers Pi, g* > 1 and let C be an open subset ofMC such that the 
functions f and gi are analytic in a neighbourhood 0 /a: = 0, a = c, v = 0, for all c E C. 
Then there exists an n-parameter family of solutions of the system such that 

ai(x) = Ci + O(xP^), bi = Oix'i^), 
where ai{x) and bi{x) are defined for c eC, |x| < a;o(c) and are analytic in x and c. 

This theorem allows us to parametrise the family of solutions near a singular point of 
a set of ordinary differential equations. We need to take each boundary point in turn and 
transform our held variables so that the held equations are in the form required by Theorem 
121 After that, it is elementary to verify the forms we have chosen for our expansions of the 
held variables near the singular points fHSl SHI [M|) . 

A. Power series results atr = 0(fc = l only) 

For the boundaries r = rh and r ^ 00 , we saw that physicality requirements hxed all of 
the boundary conditions in a relatively simple way. In the case of solitons we anticipate the 
situation will be much more intricate, as it was for the purely magnetic su{N) cases studied 
for hat spaced and for asymptotically adS spaced, and computation of lower order terms 
conhrms this assertion. This is due to the complicated intercoupling between the gauge 
functions. Therefore we now study the boundary conditions at the origin in detail, which 
are obtained by assuming power series which are good near the origin and requiring that the 
metric and held equations are regular there. Our strategy at hrst then is to take the held 
equations and convert them from diherential equations in the held variables into recurrence 
relations in the series expansion parameters. We then attempt to hnd a consistent solution 


(57) 
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to these relations. 


We begin with the held equations in the forms flTSl dSl [211 ED]), multiplied through by 
factors of /i and S as appropriate, anticipating the later series expansion. Note that we are 
here using /i in the form 


/i(r) = 1 


2 m{r] 


^2 

+ p- 


(58) 


Then we use the su(2) embeddings to hnd a rescaling of all quantities, as follows. First we 
dehne the quantities 


7i =j{N- j) (for j e {1, 2,..., iV - 1}), 


X2 — 
Ajv = 


N{N"^ - 1 ) 


(59) 


1 


and then we rescale: 


r = Aatx, 


m(r) = Xpfm{x), I = XnI-, 


^jA) = ^jAnx) (for j G {1, ...,N - 1}), (60) 


/i(r) = fi{XNx), 


S{r) = SoS{Xr,x). 


Finally, for ease of computation we dehne the quantity 


- iV - 1 + 2j. 


(61) 
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All of this brings the held equations into the forms 



- (^G + P) + i 

^2knGS (/i=s) + 


X =n 

dx^ 


dS _ 

X— -26 

dx 


X- 


d^j 

dx 





d^Uj 

dx"^ 



(~ 6 ^ 
y,n-KMXP+j'j 


duj x’^n'i] duj 
dx 2 S'q dx 


fiS'^ \%x^ ~2 

[Qj+i - Qj) Uj- 


4Si 


-Uj£j , 


(62) 


with 


4 du. 


G — kn Ij 


j=i 

^ ^da 


- 


dx 


= K,j\fG, 


h = E . , 

j=i \ dx 


J \ 

= knV, 


Kn 


N 


i=i 

Af-l 

C = E 

i=i 


(63) 


Requiring that jl, ^ and the held equations themselves are regular, and noting our rescaling 
of S, leads to the following requirements for the lower order terms: 


Sq = 1, Ujfi = 1, rho = mi = m 2 = Si = Uj^i = = 0. (64) 

(For clarity, we’ll now drop tildes.) This means the power series will have the basic forms 


00 00 


m{x) = E ^kX^, 

S{x) 

= 1 + E SkX^ 

k=3 


k=2 

00 

Uj{x) = 1 + E Uj^kX^, 

8 j{x) 

00 

= E £j,kx'"- 


k=2 k=l 


(65) 


To ‘sweep up’ factors of /i and S which came from multiplying through originally, we make 
the replacements 

fiS‘^ = l + M, = 1 + M, (66) 
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so that all the summation terms are within M and M, which are given by 

OO OO OO k — 2 

m = j2 =5^ (/ifc + 2^fc) + EE {SiSk-i + 2fiiSk-i) 

k=2 k=2 k=4 1=2 

OO k—4: k—l—2 

+ E E E 

k=Q 1=2 p=2 

OO CO CO k—2 

M = J2 ^kx’^ =^{Sk + 2fik) X^ + EE (^fiifjjk—i H“ X 

k=2 k=2 k=4: 1=2 

OO k—4: k—l—2 

+EEE fjjpfxiSk—i—pX . 

k=e 1=2 p=2 


(67) 


It can be seen that writing out the field equations using these quantities makes the picture 
considerably clearer: 


dm 

dx 


— KjV (p-G + i^) + M 


(-^ + kat (/iG + P)^ 


X^fiTl 


dS 2 knSG - ( dS 2 knSG\ ( 
— =—--h M -^--- H-7— 

dx X \ dx X J 2 Sqx’ 


( 68 ) 


9 d?Uj 


= — 2 \ m — KnXP + 


M I 

dx^ 


X 

-2 


duj 1 , , 

^ x^\ duj 

m-K„xP+-p — 


-7t (Qj+l - Qj) Uj + 


X^flTj dUj 


x^X% 

d ^2 


Uj£j , 


(and the Sj equation in fl62l) is unaltered). Thus it may be seen for each of fl68|l that the 
terms divide up into three groups - terms which have been examined previously^, those 
same terms multiplied through by a factor of M or M, and some extra terms in respect of 
the non-zero EGFs. This makes the analysis a lot more tractable. 


Now we inspect the recurrence relations given by substituting the power series fl65ll into 
the field equations fl68l) . We begin with the Einstein equations as they are much easier. 
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The recurrence relation for ruk is given by: 


[k + ( G, + + a) + ^ ^ ^ 


k-2 r 

E 

1=2 


—2mi+iGk-i — 


mi+iVk-i 

2 SS 


+M; (A; — / + l)mk-i+i + /tAT ^-h Gk-i + Pk-i 

k-l-l 

2 ^ ^ TfipGk—i—p+i 
p=3 


(69) 


where 


Vk = EEi'+ix*^' — / T 1 )q;j 
i=i 1=0 
N k-1 

g*=EE('+i)(*^ - I + l)Uj^i+iUj^ 

j=l 1=1 

N-l k-l / k-l-l 

Ck = I Sj,l^j,k-l + 2 

j=l i=l V P=1 L 

, AT k 


j.k—l—l’) 


/c—/+I 5 


fc—/—p—1 


^j-,p^j,l^j,k—l—p “ 1 “ ^ ^ ^j,r^jj,1^j,k—l—p—r 


r=2 


, (70) 


^’‘EEE 


A;—/+27 


j=l 1=2 


k-2 


Qj,k ‘^'^j^j,k ‘^'1 j— 1^j— 1 ,k T ^ ^ j,i^j,k—l 'yj— 1^j—1,1^j—1 ,k—l) 1 


1=2 


we note that Gk, Ck, Pk and qj^k are non-zero only for A; > 2 (though T]k ^ 0 for all A; > 0). 
Note that we leave rjk in terms of aj^k ~ this is purely for simplicity, as we recognise that 
Sj^k = C(j,k—0(j+i,k for all A; G N and thus terms of order P in aj depend only on terms of order 
P in £j. Upon examination, we see that ruk+i depends only on m3,..., rrik-i, So, S 2 , ■■■, Sk- 2 , 
'Uj 2, ■■■,'Uj^k and £j^i, •••, £j,k—i- 
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The recurrence relation for S is: 


kSk —2KNGk + ^KNGiSk-i + Ml [—{k — l)Sk-i 

0 1=2 ^ 


k-2 


k-l-2 

+2kn I Gk-i + SpGk-i-p 

p=2 


(71) 


and we readily observe that Sk depends only on ms, Sq, S 2 , ■■■, Sk- 2 , ^^, 2 , 

and Sjs,£j,k-i- So as long as we can solve the Yang-Mills recurrence equations to hnd a 
consistent regular solution at r = 0, then we can use those parameters to hnd mk and Sk at 
each order. 


Now we come to the much more complicated Yang-Mills equations. They can be expressed 
as the following: 


[A} - k{k + 1)5}] = z\ 

[A] - k(k + l)a‘] = b‘. 


(72) 


We dehne 5* as the Kronecker symbol, A® as the {N — 1) x [N — 1) matrix given by 


A- = (2S‘ - yy, 


(73) 


and we have dehned 2N — 2 vectors of length (N — 1) (in components): 


Si — {Si^k,S2,k, ■■■,SN-l,kY', ul — (Ul,fc, U2,fc, •••, MAf-l,fc)^; 


(74) 


and two length-(Y — 1) vectors z® and 6® representing the right hand sides, given by: 


Zi. = 


k-2 

E 

l=i L 


N-l 


+ Y, 


j=i L 


k-l-2 




—l—p—1 


p=2 


+s,+i ([y -{k-i- i)(i -21- 1)5‘] q_,_, 

k—l—p—2 


(75) 


k-l-2 




p=2 


r=2 
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and 


bUi 


k-2 

E 

1=1 

1 . 


N-l 


U 


i+iM.k,i + '^k+i^j'^k-i 

i=i 

k-l-2 


+ 2“?+i E - ^'+1 (W - (^ - ^ - 1)^}] 

p=2 


u 


k-l 


k-l-2 

-E 

p=2 

X 2 


^ ^ k—l—p—2 

— AjUpUf^—l—p-\-UpAjUk—i—p-\-—Up ^ ^ Aj^r'^k—l—p—r 

r=2 


^uiui 


k-l-2 
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+ E 


fc—/—p —1 


p=2 


(76) 


with A4,i being the {N — 1) x {N — 1) matrix given by 


■Me,/ =21(1 + l)mk-i+i + ^ ((fc — 21 — 3)Sk-i-2 ~ (k — l)(fc — ‘^)^k-2) 

k-l-2 

2 ^ ^ l(p I 1) *S*p'n7.^_;_p_|_i, 

p=2 


(77) 


and being the (A^ — 1) x (N — 1) matrix given by 


■Aik,I =(l + 1 ) 


2(1 - l)mk-i - E<5i_2 + 2KNPk-i-i 


1 /, /!> I Vk-l-3 Vk-l-5 

+ 2)Mk-l-3 + 25-2 + 25 ' 2£2 


(78) 


k-l-2 

E 

p=2 >- 


2(1 \')MpTnk—i—p 2tviq MpPk—i—p—\ 


^p+ihfc—/—p —3 

45^ 


The system with which we are left thns bears similarities to the pnrely magnetic system^, 
and so a sensible strategy here is to nse resnlts in that work. We expand the matrix 
and the solntion vectors u* and £* in the left- and right-eigenvectors of A*, and we attempt 
to prove that each new vector in the expansions reqnires only one new parameter each, 
and hence that the solntions to this eqnation are consistent in their choice of parameters. 
Therefore we are in a position to state onr proposition: 


Proposition 3 The system given by the recurrence relations fl69p . fl7ip and the mutual 
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Yang-Mills system 


{A - k{k + l)I)uk+i = bk+i, 
{A - k{k + l)I)Sk = Zk, 


with conditions 


cr^Sk = dkOki 
cr'^ut+i = 4 / 3 s, 


(80) 


has a consistent solution in its parameters, i.e. that all parameters mk, Sk, Uj^k, ^j,k depend 
only on m 3 , ...,mk-i, Sq, Sk-i, Mj, 2 , ■■■,Uj,k-i and 8j^i, ...,£j^k-ii thus each new parameter 
in the expansions can be obtained from those previously calculated. That is to say, each 
order of the gauge field expansion will contain only one new parameter, for 1 < k < N — 1; 
and since the metric expansions are determined entirely by physicality and the gauge field 
parameters, this gives us therefore 2N — 2 solution parameters in total. 


Proof We have already stated that if the gauge functions have a consistent solution then 
the metric functions will too, so we focus on the Yang-Mills system fl7^ ISO]) . The key is 
that because both systems contain the matrix At, we may appeal to lemmata 1, 2 and 3 
formulated by Kiinzle^, the main results of which we summarise below: 

• The matrix At has eigenvalues 1.2, 2.3, ... , {N — l).iV; i.e. the eigenvalues of At are 
k{k + l) for k e { 1 ,..., Y- 1 }. 

• The left and right eigenvectors of At, {cr^} and {vk} respectively, are given by 

= wziiA^ vi = ^QkU - 1), (81) 

where Qj{r) is dehned using Hahn polynomials, with normalisation 

N-l 

(^ 2 ) 

i=i 
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and normalisation constants dk given by 


(Ar +A;)!(iV-fc-l)! 


(83) 


(N - l)!(Af - 2)]k{k + l){2k + 1)' 


(Fnrther details on this may be found in the literature^.) 


• The system given by 


{A - k{k + l)l)uk+i = bk+i, 

= 4A, 


(84) 


has a consistent solution near x = 0 regular in all held variables. 

Examining the Sj system in this notation, we hnd it is very similar to the form of the Uj 
system: 

{A-k{k + l)l)£k = Zk, 
cr^Sk = dkOk- 

Therefore lemmata 1 to 3 will apply equally well to the mutual system flSTl [85]) - for example, 
see the proof of local existence near r —)■ cxo due to Kiinzle^ - and so we can follow the basic 
form of this proof. 

Examining the right-hand sides of fl75l 1761) and hnd that 


bi = b 2 = zi = Z 2 = 0 . 


( 86 ) 


Now, examining lemma 3 in Kiinzle^, we hnd 


U2 = AUi, U3 = ^2^2, 
S\ = 9iVi, 82 = d2V2'-, 


(87) 


and by a similar argument, we get 



( 88 ) 


for 3 < fc < — 1. Here, 8 ^ 


k 


and are special solutions hxed by the requirements 


cr^ ■ bk+i = ■ Zk = 0. 


(89) 
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We emphasise that the £]^, u]^ are entirely specihed by using the results in Kiinzle^, 
and then the 2N — 2 arbitrary constants f3k, dk entirely specify 81 and 


To perform the analysis it is convenient to expand our equations in terms of the left and 
right eigenvector basis, and we may use the same basis for both sets of equations according 
to the above facts. In this basis A* may be written 


N-l 




a=l 


and the gauge held vectors can be expanded in the same basis as 


(90) 


N-l 


N-l 


'“Ui = E ulvi si=T. m, 


1=1 


1=1 


(91) 


for constants Ul and E^. Note also that 


and (ITID imply that 


U^k=f3k. El = ek. 


(92) 


Finally, for later convenience we dehne the quantities 


N-l 

drs = da^ E 

i=l 


N-l 


dr St = da^ H <vivivi. 


2=1 


We may thus express fl72|) as 


(93) 


k-2 


k){a + k + l)El — — 


N-l 


Nk^iE^ + 2a(a + 1) El_,_^ 


1=1 L r,s=l 

k-l-2 N-l 

+ E “(“+1) E 

p=2 r,s,i,u=l 

+^/+i ([a(a + l)-{k-l-l){k-2l- 1)] 

k-l-2 /N-l 

+ E E 2“(o + 

p=2 \r,s=l 

k-l-p-2 N-l '' 

+ E E 

r=2 q,s,t,u=l 


( 94 ) 
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k-2 


k){a + k + l)U'i^ = -Y,[Mk,iU^ 

1=1 

N-l . . 

+ E + 1 ) + + 1 ) 

r,s=l ' 


k-l-2 N-l 


+2 E E (“(“ + ^Kudr.u;u;_,ui-i-,-i) 

p=2 r,s,t,ii=l 

-Mi+i [{a-k + l-l){a + k- 


k-l-2 /N-l 


p=2 \r,s=l 


E E +1)++1) dt.u;-NL,-,-i 


k-l-p-2 N-l \ 

+5 E E iv{v + i)d:Au;_N‘r-iU^_,_,_,_,)\ 

r=2 s,t,u,v=l / 


.2 A^-1 / N-lk-l-2 

+1# E e\eu_a + E E ^Ueie. 

0 r=l V s=l p=2 

Now, for the lower order terms we note from fl9^ that 


k—l—p—l^r^s 


(95) 


£! = 9i, B," = 92, C5;=ft, U^2=I32, 


(96) 


and from the right hand sides of 


we see that 


= 0 for a > 1, = f/“ = 0 for a > 2; 


(97) 


SO that the hrst two terms in each gauge held expansion are dehned by one parameter each. 
What we wish to prove is that this is true in general, and that every new term in the 
expansion is determined by one new parameter each, until N — 1 parameters have been 
introduced for each gauge held (note that this is not assumed - it emerges naturally from 
the size of matrix (173|l ). Therefore we prove the following: 


Proposition 4 = 0 for a > k. 


Proof In order to show this, we use a form of 2-proposition ‘ladder’ induction, dehned by 
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the following sentences: 


E{k): 


= 0, 

for a > u, for each z/ = 2,., 

..,/c 

U{k): 


= 0, 

for a > u, for each z/ = 2,., 

..,/c 


( 98 ) 


Our strategy is as follows. 

• We know that E{2) A U{2) is true. 

• We make the induction hypothesis E{k) A U{k) for some kgN, 2<K<iV — 1. 

• We show that E{k) A U{k) E{k + 1). 

• We show that E{k + 1) A U{k) U{k + 1). 

• Thus we have shown that E{k) A U{k) E{k + 1) A [/(k + 1) and the induction is 
complete: E{k) A U{k) is true for all k >2. 

We will need to make use of lemma 4 in Kiinzle^, that is, 


= 0 if r + s < a, (99) 

in order to endure that strategic terms of flMj) and fl95H vanish. In addition to this fact we 
will need another analogous fact: 

d-rst = - ^)^Carst = 0 if r + s + t<a (100) 

where Cijki is dehned below. We shall now briefly prove this fact, by way of a lemma. 
Lemma 5 

7V-1 

Cijki = - f)Qj{r - l)Qk{r - l)Qi{r - 1) = 0, (101) 

r=l ' ' 

if i + j + k < 1. 


Proof First we note that Cijki is symmetric on all of its indices. We shall use the following 
induction argument: 





• We show that Cuu = 0 for 2 + fc < /; 

• We make the induction hypothesis that 1 + j + k < I implies cijki = 0, for all j G 

{1,J} for some JgM, — 1; 

• We show that Cuki = 0forl + J + A;<Z implies that Ci^j+i^ki = 0 for 2 + J + /c < Z; 

• Finally we note that because Cijki is symmetric in its hrst two indices, we are done; 

since the line above implies that if = 0 for / + 1 + Zc < Z, then = 0 for 

/ + 2 + Z^ < Z, and therefore induction over one index is induction over both indices. 
Hence we conclude that 

Cijki = 0 if Z + j + Zc < Z. (102) 


Firstly, note^ that Qi{r — 1) = 


N-l 


(103) 


e.g. if 2 + Zc < Z. 
Next, assume that 


N-l 


(104) 


if 1 + J' + Zc < Z, for allj G 1 < J <A^ — 1. Now we use the recurrence relation for 

the Hahn polynomials 


- rQi(r) = diQi-i{r) + {k + di)Qi{r) + hiQi+i{r) (105) 


where 

(* + l)(iV-z-l) (,-i)(iV + z) 

2(2z + l) ’ * 2(2*+ 1) 


(106) 


in (11041) : hrst with Qi{r — 1), then with Qk{r — 1), then we let j i—)■ j — 1 and equate the 
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expressions. This yields 



h-i 

+ {bj-i + dj-i — bi — + biCij-i^k,i+i] ■ 


(107) 


Now to complete the induction step, we let j = J + 1, with the restriction 2 + J + k < I, i.e. 
J + k < I — 2, and consider each term of (11071) : 

• = 0, since 1 + (J + 1) + A: < /, i.e. J + k < I — 2; 

• = 0, since 1 + J + fc < /, i.e. J + fc < / — 1, and we know J + k<l — 2<l —1] 

• = 0, since 1 + J + k < I, i.e. J + fc < / — 1, and we know J + A:</ — 2</ — 1; 

• ci,j,k,i+i = 0, since l + J + k<l + l, i.e. J + k < I, and we know J + k<l — 2<l. 

Therefore, all terms on the right hand side of fllU7p vanish and hence Cijki = 0 G 
{1,N — 1} if 1 + j + k < 1. Finally, using the symmetry argument described earlier, we 
can say that Cuki = 0 ^ iii + l + k < I, and hence the lemma is proven. □ 

Corollary 6 This means, using (ISTl) . that we can write 


'Kst ^0 if T + s 1 < a. 


(108) 


Thus we have proven (11001) true. 

So now we are equipped to take fIMl) and in turn, making the appropriate induction 
hypothesis E{k) A U{k) for some 2 < k < iV — 1; hence we now substitute k = k + 1 into 


each term of fIMl) and prove that E{n) A U{k) E{k + 1). 


• ^ N'k+i,iE^ = 0 since 1<k — 1<k + 1. 


i=i 


• 2a(a + 1) E E E d^sU[E^^_i: 


1=1 r=l c=l 


Now r + s<K<K +1, SO that d^^ = 0 here. 
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K—Ik—/—I p—1 I K—l—p 

•<*(« +1) E E EE E 

l=l p=l r=l s=l t=l 

Here, r + s + t<K — 1 <k + 1 so = 0. 

K—1 

• >S'i+i(a + k — l){a — k + I + = 0 since k — I < k — 1. 

1=1 

K—l K —1—1 p—1 K — l—p 

. 2a(a +1) E s,+, E E E d‘,u;_,EU-y 

1=1 p=2 r=l s=l 

This is similar to term 2: we get r + s<K — I — 1 <k — 2, hence 6?“^ = 0 here. 

K—1 K—l—lK—l—p—lr—lp—lK—l—p—r 

. a{a +1) E s,+i E E E E E 

1=1 p=2 r=2 q=l s=l t=l 

Now q + s + t<K — I — 2 <k — 3. Therefore = 0 here too. 

Hence we can say that E{k) A U{k) + 1). Now we must take (j^ to prove the 

other statement: that E{k + 1) A U{k) U{k + 1). 

It is clear from (j95ll that the left hand side and the first two terms on the right hand side 
are identical to those previously examined^, so these are already proven to be zero under 
this process, especially as that examination^ makes no mention on the form of A^y/c and 
so the appearance of i is unimportant^. The next three terms are very similar, but we’ll 
examine them anyway. 

K,— l 

• — k + I + l)(a + k — l)U^_i = 0 since k — I < k — 1. 

1=1 

• E A+i “e ' E ’‘e’’ (h(o +1) + 0(0 +1)) d‘rKvu-r 

1=1 p=2 r=l s=l 

Now r + s<«: — / — !<«: — 2, SO (i“g = 0 here. 

K—1 K—1 K—/—p—1 p—1 r—1 K—i—p—r W—1 

• 5 E E E EE E E v{v + 

1=1 r=2 s=l t=l u=l v=l 

Here u + t < k — I — p — 1, therefore = 0 unless v < k — I — p — 1. But then, 
v + s<K — I — 2 <k — 3. So dg^ = 0. 


The final two terms are less familiar but present little problem. 


X2 


K-l N-l 


V- 


-[•dr 


In this term, note that we have E^_i in the product. But I > 1 and so k — 1<k — 1 
and this term is zero. 
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A 


N 


E E E 


452 — 

'J l=l r,s=l p=2 


Finally, we note the Up_i multiplicand in this term. However we also have p — 1 < 
K — I — 2<k — 3, so once again this term is zero. 


Hence we have shown that + A 17 (k) U{k, + 1). Thus the induction is complete, 

and Proposition m is proven. □ 

Therefore, we can now dually express Uj and Sj as dnite power series, in vector form: 


N-l 


u{x) =Uo + 

k=l 

N-l 

= E &k{x)x^, 


k+1 


(109) 


k=l 


where Uq = (1,1,..., 1)^ is an — 1-vector and 


N-l 

I3k{x) = f3k{x)vk = Y. UfvkX^-’^, 

l=k 

N-l 

Gk{x) = 6k{x)vk = E^VkX^-^. 

l=k 


( 110 ) 


We dedne 2(iV — 1) initial parameters {/3k, 9k} where 


^k = /3k{0), 9k = 9k{0), (111) 

in which the constants and Ei depend quite complicatedly on the constants {^k,9k}- To 
obtain the boundary conditions fITSD then, we substitute Uj back into the expression for uij, 
and augment them both with possible higher order terms. Therefore we have shown that 
the gauge held equations are entirely specided by 2N — 2 initial parameters which can be 
chosen in a consistent way. 

Finally, we dnish od our proof of Proposition [3] with the following argument. Assume 
that for some k we know all the parameters m^, ...,mk-i, So,Sk-i, Uj^ 2 , ■■■,Uj,k-i and 
Sj^i ,..., £j,k-i- Using the recurrence relation fl6^ we can certainly work out m^. The detailed 
study of the Yang-Mills equations in the form 0721) has shown that we can work out Uj^k and 
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Sj^k- Finally, we can use the final recurrence relation fl7T|) to work out S'a,, and Proposition [3] 
is proven. □ 

Before we go on to prove local existence near the origin, we note an interesting consequence 
of flllOp . It may be noted that the dependence of both j3k{x) and 6k{x) on x is the factor 
x^~\ hence we might expect that in fact I3k{x) oc 6k{x) for each k. This is fairly easy to see 
once it is realised that we have expanded both gauge fields in the same basis of eigenvectors, 
but to be precise, we can prove a minor lemma that will help confirm our results later. 

Lemma 7 

Uk+i = TkSk for N — 1 constants r^, (112) 

for k G {1, ...,N — 1}, and hence 

f^k = nOk, (113) 

no sum over k. 


Proof Take the second equation in each of flH^ and fl85l) : 

cr^Uk+i = dkl3k, a-^Sk = dkOk, (114) 

(for dk defined in ([83])) and define Tk = ^ (no sum on k). Then it is clear that 

cr’^Uk+i = TkfT^Sk. (115) 

Since this is a scalar product and we have not mentioned the form of at all, this must 
hold for each pair of vectors Uk+\ and Sk, i.e. 


and immediately. 


(no sum on A;).n 


— XkSki 


Pk = 'TkOk, with Tk 


fdk 

Ok 


(116) 


(117) 
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Notice this is not saying that the functions Uj and Sj are multiples of each other; rather 
that at each order, the expansion vectors will be at most a scalar multiple different from 
each other. Later, this gives us a form to aim for when rearrange the equations to make use 
of Theorem [21 


B. Local existence of solutions at the origin r = 0 {k = 1 only) 

Now we can use the results obtained in Section m to prove the existence of solutions 
regular at the origin and analytic in their initial parameters. First we must rewrite the 
equations in a form applicable to Theorem |2l Hence we will state our proposition: 


Proposition 8 There exists a 2N — 2-parameter family of solutions to the field equations in 
the form / fldl \TU. [^]) which are regular at r = 0. These solutions are analytic in some 
neighbourhood of r = 0 and in their initial conditions {/3fe,0fc} fillip . 


Proof We begin with the scaled equations from Subsection IIII Al fl62p {sans tildes), except 
that we take the original dehnition of S so that the hrst term of the expansion is Sq. 

It is best to begin by examining the equation for the MGFs. Introducing the vectors 
u = (mi, M 2 , ..., mat-i)^, £ = {£ 1 , 82 , , we may write this equation as 


2 d‘^u 

X li-rr + 2 
dx^ 


^3 ^3 


+ Iw + llLv = 0- 

dx 2 4^52 "• 


(118) 


introducing two N — 1 vectors 


W =(>Vi, >V 2 ,..., WN-if with >V, = (g,+i - qj) Uj, 
V =(Vi, V 2 ,..., V^_i)'^ with V, = u,£]. 


(119) 


We notice that this is very similar to the corresponding purely magnetic fornt^, so parts of 
the analysis will carry over similarly. Using fITSp . we can write kV, as 

Af-l 

Wj = 2uj — ^ UjAtul (No sum on j). (120) 

i=l 
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Also, we make use of the fact that as functions dehned with the eigenvectors of A*, 
and Ok{x), will satisfy 

Af3,^{x) = k{k + l)l3f.{x), A 6 k{x) = k{k + l) 6 k{x). ( 121 ) 


We should state that we begin by assuming nothing about the properties of the functions 
f3i^{x) or Ok{x) at this stage; for instance we do not yet assume they are regular at a: = 0 
nor that they are of the form (IllOp . 

Multiplying by the left eigenvectors cr^ (which we assume to be constant and independent 
of x) gives: 


dx^ 


+ 2 


a; X ri 
m-xnNP + -- — 


A. 

dx 


[cr^u) + -cr^W + 


A 


N 


X 


4/iA2 


cr^V = 0. (122) 


Using the orthogonality of eigenvectors, whence a^(i^{x) = 0 iff A 7 ^ 5, we dehne — 1 
functions cr^(3i^{x) = ^l{x) (no sum on k). Using (1109^ . this transforms (11221) into 

0 =a:V ( + 2 {k + l)x’^^ + k{k + l)x'^-^^l{x)) 

\ dx dx J 

+ 2 (^m - Kf,xP + ^ - + {k + l)x'‘(l(x)j (123) 

+ + 1^,7‘v. 

2 

The hnal term involving V does not cause concern since it is of quite high order; the compli¬ 
cated part of the analysis is resetting W in a suitable form. Fortunately we may appeal to 
two lemmata courtesy of Kiinzle and Oliynyk— (see equations (7.3) and just below (7.7)): 

Lemma 9 

.. N-l z 

-W = -^fc(/c + l)/3fc(a;)a;^+^ + ^y;x' (124) 

k=l 1=2 

for vectors Qi and some integer z > 2. □ 

The details here are unimportant (we shall show an analogous proof later), as are the 
precise value of ^ or the forms of grp the key point in this result is that W ~ O(x^) at least. 
In addition, we have 
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Lemma 10 

(T^Ql = 0 i/ / < fc + 2. □ 


(125) 


Concerning the final term in fll23p . we may establish that 


N-l N-1 


^ 


.X 


l-\-p—2 


(126) 


1=1 p=l 


If we introduce new variables 




dx 




(127) 


divide (I123p through by rearrange and use (112411125p we can finally rewrite (I123p as 


= -2{k + l)ipl - k{k + l)^lx ^ + k{k + 1)^^ 


X 


-1 




2 / r> 

m - KmxP + — 
X^fl \ P 


x^rj 


4^ 


AT-l N-l 

ii E E 

1=1 p=i 


X 


4521 ((*= + XI + ^i’l) 

i+P-2 _ ^-1 ^ 


(128) 


<r'‘gi+k+X- 


l=D 


Now we leave the Yang-Mills equations for a while as we must consider the expansions 
of the various quantities G, qj, P, rj and (. In order to do this we must note that the first 
eigenvector Vi = (1,1, This is easy to prove if we observe that 


Y - 1 

v\ = — -: 2 -Fi(—1, —i + 1; —N -f 1; 1), (129) 

IS — % 


and then use the special case of Gauss’ theorem; that is, 

2 -Fi(— n, b; c; 1) = ^ for 91e(c) < 91e(—n + b), (130) 

(c)(n) 

for n G and where (a)(„) is the Pochhammer symbol defined by (a)(o) = 1, (a)(n) = 

n—1 

n (a + i) in> 0 )^. This means that all Uj ^2 are equal for all j, and similar for Hence 

i=0 

we write 

Uj ,2 = /3i{x), £j^i = 9i{x), (131) 
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(for all j) and we introduce the expansions 


/3i{x) = /3o + xPi{x), 6*1 (x) = 9o + x9i{x), 


(132) 


where again, we assume nothing about the functions /3i(x) or 6^1 (x); they can be written 
in terms of and but the precise expressions are not important. Note also that from 
now on, a hat" indicates a function we have introduced to deal with higher order behaviour, 
which is analytic and regular in all transformed held variables with which they are dehned. 
Thus we obtain 


G 

P 

c 

V 


4/^0 I Ji, 


Kn 


W + x"G, 


2{N - 1 + 2j)x^/5o + x^qj, 


^x^ + x^p, 
Kn 
02 

^x'+X^C, 

Kn 

2k 


+ XT). 


N 


(133) 


Now we can hnally tackle the Einstein equations in fl 68 ll . Using (I133p . they may be written 


dm 


3 



8kn J 


dS _ 

X— = X Da- 

dx 


S , 4 ^ 

X + X m, 


(134) 


The equation for S is now in the required form, and we note by the theorem that this implies 
the expansion is of the form S{x) = So + O(x^), a fact we will need momentarily. 


If we dehne a new variable for m by 


1 9"^ 

X = — (m- msx^) with m 3 = 2/3o + —2-, 
X oAv 


(135) 


then we hnd 


dx 

X— = —3x + xm. 
dx 


(136) 
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Note also that /i = 1 + {21713 ~ 2% + ^). This equation is now also in our desired form, 

and from Theorem [2] we get the behaviour x = 0{x) and hence m = msx^ + 0{x^). 


Now we can go back to the MGF equation fll28p . Using fll34p . and writing /i ^ l + 


S ^ = S'o ^ + x'^Sb, and S = Sq^ + x'^Sc, we get 


-2 _ c -2 I ^21 


tjC^ tjC^ T1 \ 

-2 ( m — xknP + — — TT^ ] = and 


£2 4^2 




N-lN-1 




(137) 


^ (=1 p=l 


Therefore we can write fll28p as 


= -2{k + l)V’fc + xH2,k - cr^«/fc+2- 


dx 


(138) 


We make one hnal change of variables here, as we let 


and hnally get our equations in the desired form: 

= -2(/c + l)'4)l + x'H2,k- 
dx 


(139) 


(140) 


Now we turn our attention to the equation for £j in fl6^ . We shall need to dehne analogous 
transforms to those for and ipl, that is 

(T^0k{x) = ^l{x) (no sum on k), (141) 

We recall that lemma [7] established that = TkOk for some constants r^. Therefore, we 
can deduce that = Tkil i’l — Xk'ipl'^ and hence if and ^ are regular and analytic 

near a: = 0, then so should and ^ be, and they should exhibit the same behaviour at 
the origin as and as dictated by Theorem [2l We do not feel that completes the proof; 
rather it gives us a method of conhrming our result for the analogue of (I140p for the EGFs, 
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since it means the equations for the two gauge helds should look almost identical. 
We begin with the EGF equations in the following vector form: 




d‘^£ 

dx"^ 


= XfJ, 


X dS 
S dx 


- 2 


dS 

dx 


+ -Z, 


(143) 


with Z in component form as 


Zj — 2'yjUjSj 


(144) 


We proceed as before: we multiply through by the left eigenvectors cr^ and substitute in our 
new variable (1141 p . which gives us 


0 


=a;V {+ 2kx^-^^ + k{k - 

\ dx^ dx J 

\b dx J \ dx J 


(145) 


Similar to before, the troublesome part is writing .Z in a form that renders it applicable to 
Theorem [2j but similarly, it may be expressed in terms of the matrix A* fITsp . and thus we 
may formulate two lemmata analogous to [9] and [10], found in Kiinzle and Oliynyk— : 


Lemma 11 


N-l y 

k=l 1=3 


for some constant y G N"*", y > 3. 


(146) 


Proof We may write Z in components as 

N-l 

zj = ( 147 ) 

i=l 

We let Mj = Mj — 1, so that the lowest order term in Ui is order x^. Thus 

N-l 

Zj = (148) 

i=l 
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Consider 


N-l 


N-1 N-l 


E = E E 

i=l i=l k=l 

Using the second equation in fll2ip . we obtain 


(149) 


N-l 


N-l 


E = E 


(150) 


2=1 


k=l 


Finally we note that the other terms in fll48p are at least of order x^, and therefore our result 
follows. □ 

Once again the value of y and the form of the vectors fi are unimportant compared to 
the form of fll46p . We also prove the following: 

Lemma 12 

a^fi = 0tfl<k + 2. (151) 


Proof Following the previous lemma we may consider Z fll48p as a sum of three parts: 


N-l 


N-l 


N-l 


Zj =+2 E +E 


(152) 


2 = 1 


2 = 1 


2 = 1 


We use the previous eigenvector expansions of the matrix Aj fITsp . the dehnitions 
IIOOP and the second equation in fll2ip to give (in components) 


ajZj =k{k + l)^l{x)x 


2N-1 1-2 


+ 2k{k + l) ^ '^dkdpi_p_-^ep{x)/3i_p_i{x)x^ 

1=2 p=l 
3V-1 i-4 l-p-2 

+ k{k + l) EEE dkd!l^p^l_p_q_2eq{x)Pp{x)Pl_p_q_2{x)x^ 

1=3 p=l q=l 


(153) 


(where we note that ^k{x) and 6k{x) are zero for k > N). Comparing this equation to fll46p . 
it is clear that we wish to establish that the last two terms are zero if Z < fc + 2. Consider 
them in turn. The second term contains d^ i_p_i, which by fl99p is zero iip + l—p— l<k, 
which it is if Z < Zc + 2. The third term contains cZ^p which by fllOOp equals zero as 
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long asq + p + l—p — q — 2<k, i.e. I < k + 2, again implied by / < A: + 2.n 


Therefore we may write: 

= k{k + 1)4V + ^ (154) 

l=k+2 


Similarly to before, we substitute in tfjl fll42p . divide through by rearrange fll45p and 

use fll54p to get 




X 


-1 


p 


y-k-l 


p 


1 dS 


■' ^'"fl+k+l^' + -g-^{x'fk^Hl) ■ 


l=l 


(156) 


Again letting p ^ = 1 + x'^p, using S ^ = S'q ^ + x'^Sc, and substituting the expansion for 
the second Einstein equation fll34p . we obtain 


X 


dx 


—2{k + + xSk- 


(156) 


This is now in the form we require. As a last step we must rewrite fll27p using the new vari¬ 
ables fll39p : this transformation does not signihcantly affect the structure of the equations, 
which become 


Hence, all our held equations are in the required form we need for them to apply to Theorem 
m and furthermore, it can be seen that the result fll56p is indeed almost identical to the 
other transformed gauge held equation fll40p . just as we expected from Lemma [71 
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So, to summarise, we have brought the field equations into the following forms: 



( 158 ) 




It may be verified by substituting back that the functions m, Sa, Gk, 'H2,k and J'k are all 


regular in all transformed field variables at x = 0. At last, using Theorem [2] gives us the 
following expansions at the origin: 


X(x)=0(x), ^l{x) = ^Iq + 0{x), ijl{x) = 0{x), 


(159) 


^(x) = ^0 + O(x2), ^l{x) = ^lo + 0{x), ijl{x) = 0{x). 


By transforming our variables back, it can be checked that this gives the required behaviour 
for each of the field variables near the origin, completing the proof of Proposition (Hid 

C. Local existence of solutions at the event horizon 

Now we turn our attention to the field equations in the black hole case, i.e. at r = (for 
7 ^ 0). We assume the existence of a non-degenerate event horizon, so that /i(r/j) = 0 but 
> 0 is finite. Fortunately the black hole case is a great deal simpler than the previous 
soliton case. We begin with our proposition: 

Proposition 13 There exists an 2N — 2-parameter family of local solutions of the field 
equations 051 UR l21[ f^) near r = r^ analytic in r^, A, Uj^h, and p = r — rh such that 


T{rh + p) = p'hP + 0{p^), 

S{rn + p) = Sh + 0{p), 
ujjirh + p) = + 0{p), 

aj{rh + p) = alf^p + 0{p^), 


( 160 ) 


where p'f^ and Sh are fixed by physicality requirements. 
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Proof Let us define a new independent variable^'^'^ x = r — Vh, and also define some new 
dependent variables: 

p{x) = r, 


Ai(a;) = 



A 2 (a;) = S{r), 


'ipijix) 


«j(^) 

X 




(161) 


‘■■"'-‘rn'i 


, ^2j(x) = 


The field equations take the form: 


ju(r) dujj 
X dr 


dp 

X— = X, 

dx 


dA 

X— = 0, 
dx 


X- 


dXi 

dx 


— ~Ai + J-h + xQi^ 


hi 


X 


d'lpij _ iijX2 


dx 


- ~i’ij + 


p 


2 ’ 


dX2 ^ 

= xG2,hi 


dx 


X- 


d^ 

dx 


(162) 


X 


^2ji 


X- 


__ 

dx A 1 A 2 


^ /,2 rj, j, ^ 


^2j (V’li - V’lj + l) , = -^2j + Vj^h + 


where 


(V't - di-i -k(N+i- 2 j)r-±^ 


P A2 2p3 ^ 
.7=1 


X 

N 


p 


2A 1 

E E W'l; - Xhi) 


Wm = 


J = 1 


Tj,K = - Pipy (^k-ip% + l (dl+i + Ai-i) ) ■ 


i=i 


( 163 ) 


43 


















To finally put the equations fll62p the required form, we let 


7 / 


) -^1 •^l ^2j ^2j 'Pj,h- 


( 164 ) 


This makes the non-conforming equations take the form: 

= -^1 + xQi^h, = ~'^^3 + = -^2j + xiij^h, (165) 


where 


Gl,h —Gi,h — 




dJ^h _ _ 'GiMij - i’w) dT, 

dp Ai d'G2j 
dVj,h dVj,h 
dp di)2j ’ 


A 1 A 2 




•^j,h — 2 ^2 


ilj ^ "Glj / , ,, ^ , 2 A 26 j 




Aip- 


:{'Gij ~ 'Gip+i) + 




(166) 


Hence, examining (116211 and fll65p . we can use Theorem [2] to show that there exist solutions 
to the equations (lT5l fTHl [211ESP of the form 


p = rh + 0{x), 


Ai = 0{x), 3pij = V’lyft + 0{x), = 0{x), 


(167) 


A2 = A2,o + 0 (a;), 'ip2j ='G2j,h + 0 {x), ^2j = 0 {x), 

with p, Ai, A 2 , V’ 2 i, and ^ 2 j all analytic in x, Vh, ojj{rh), ot'j{rh), A, and S{rh). 

Transforming back to our original variables gives us the correct behaviour and analyticity. 
When we £x and A and choose Sh such that S' ^ 1 as r ^ 00 , this gives the expected 
2N — 2 parameters. □ 


We have thus proven existence of solutions to the field equations for a black hole in some 
neighbourhood of the event horizon r = Vh, satisfying the boundary conditions fj48ll . 
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D. Local existence of solutions at infinity 


Now we prove existence locally as r ^ oo, which applies to both black hole and soliton 
solutions. We note that, as in the adS spherically symmetric case, it is relatively easy to 
prove existence here and we need only go to first order in the field variables, unlike in the 
asymptotically fiat case where higher order terms were needed and the analysis was much 
more involved^*^. 

Proposition 14 There exists an 2N-parameter family of local solutions of the field equations 
near r = oc, analytic in A, ojj^oo, M and r ^ such that 



(168) 


Proof We transform our independent variable to x = r and introduce new variables^i^*^ 


Ai(x) = 2m(r), = aj{r), 


(169) 


X,{x) = S{r), = n;,(r), f,,{x) = 


Then the field equations take the form: 




( 170 ) 
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where 


J^oo = -1- kx^ — x^Ai = 

3 

I N ^ N-l 

Qoo,l = “ Tf ^Ij ~ OTT \2 '^‘2j ('^li “ V’l J + l) ) 


2 ,=i 


00-2 


N-l ^ Af 


- 4- - - - ^ 2 ,+! -kiN + 1 - 2j)Y , 


i=i j=i 

N-l 


N-l 


Qoo,2 


i’ljii’ij - V’i,i+i) - 2A2 J2 4 ’ 

00-^ i=i j=i 


2 TIX 2 


^00,li - V’lj+l) - ^ 2 ,i-l(V’lj-l - ^li)) + 

c 




jf^oo,2 


Aq 


'Hoo,2i + 1 + ^2 j ^ 


cxd'^'2 




^ ^^2j(3xAi 2fc) ^ X ^2j^oo,l ^ ^2jQoo,2 




F 

n< 


A 9 


(171) 


Since l//i is at least of order as x —)■ 0, it can be observed that all of these polynomials 
are non-singular as x 0. Therefore using Theorem |2l we have solutions to these equations 
with the following asymptotic behaviour: 


Ai(x) — Ai^o + 0 ( 2 ^)! i^ij — i’ljfl + 0{x), iij — iijfi +0{x), ( 172 ) 

A2(2^) = A2,0 + 0 {x'^), lp2j = t/’2y0 + 0 {x), ^2j = ^2j,0 + 0 {x). 

Therefore we have proven local existence of solutions at inhnity, and Theorem [2] confirms 
that the functions exhibit the required behaviour near inhnity fl5T|) . Also, by rescaling the 
time co-ordinate in the metric we can £x Soo = 1 so that the space-time is asymptotically 
topological adS - this satishes the boundary conditions flST)) . and the held variables are thus 
analytic in M, x, Uj^oo, C(j,oo, Cj, dj, and A. Finally, hxing Vh and A, we get the 4A^ — 3 
parameter family of solutions that we were expecting. □ 


IV. GLOBAL ASPECTS OF NON-TRIVIAL SOLUTIONS 

Having proven local existence at the various boundaries r = 0, r = and r 00 , 
we now turn our attention to proving that those solutions may be patched together into 


46 










global solutions, i.e. solutions which begin aX r = rh for black holes (r = 0 for solitons) and 
remain regular throughout the range Vh < r < oo for black holes (0 < r < cxo for solitons). 
First though we prove a more minor proposition concerning £j^ which it is worth including 
nonetheless. 


Proposition 15 £j{r) is monotonic for all j G — 1} - i.e., £j{r)£j{r) > 0 Vj G 

{1,..., — 1} and Vr such that < r < oo for black hole solutions, 0 < r < oo for soliton 

solutions. 


Proof Inspired by Section m we may write the Yang-Mills equation fl29|) in terms of the 
matrix A* fl73|l as 


N-l / n 


N-l 




1=1 ^ ' i=l 

Post-multiplying by (for some integer 1 < C < Y — 1), using flllOl I12ip and summing 
over j gives 


N —1 / 2 \ ^ ^^ 

E ( "Y E (174) 

i=l d / j=l 

It is not important what the value of C'(C' -f 1) is - we merely note that it is positive (as an 
eigenvalue of A*) and hence we define C = C{C + 1) >0. Now what we wish to prove is 
that the above relationship holds for each summand in i, and hence that we can diagonalise 
the system, decoupling the equations. We rewrite fll74p to clarify the situation: 



Cu% 


vh = 0 


(175) 


where 0 = (0, 0,..., 0)^, the zero vector of length Y — 1. We now consider the system 
of vectors as being labelled by i and indexed by C\ i.e., we take the set of vectors 

}. Elementary linear algebra tells us that if we can prove that {u^} 
forms a linearly independent set of Y — 1 vectors of length Y — 1, it is a basis, and hence 
fll75p shows that the eauation must hold for each value of i. Considerine that the vh are 
defined using Hahn polynomials it should be no surprise that this holds in our case, due to 
inherited orthogonality relations and so on; but we shall show the direct proof in any case. 
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In other words, we would prove that given — 1 coefficients Ui, 

N-l 

^ UiVc = 0 i^i = 0 Vi, (176) 

i=l 

where for now we assume nothing about the coefficients z/j. It is simple to show that z/j = 0 Vi 

N-l 

is a sufficient condition for ^i'^c = 0- To show that it is necessary, we need to derive 

i=l 

an analogous orthogonality relationship to fl82|) . Comparing the results flHTl [821 [83]) from 
Kiinzle^ (see lemma 2) with the orthogonality and dual orthogonality relationships for the 
Hahn polynomials^, which are 


N-l 

^ ^ Qn(^)Qm(2^)p(^) 

x=0 

N-l 

^ ^ Qn{^^^Qn{y')'^n 

n=0 



1 

p(x) 


5 


X 


we hnd that in our case. 


7 r„ = d-iiV(iV-l), 


p{x) 


X -\-l 

N{N -x-iy 


(177) 


(178) 


Using these with fll77p and the above results from Kiinzle it is possible to show that 


N-l 

= ^k- (179) 

c=i 

The analogy to flH^ is obvious. So, to prove (1176!) . we multiply the left-hand side through 
by (for some 1 < < iV — 1), sum over C, and use fll79p to show that 


AZ-l Af-1 


Y ^idc^kVh = 0 

C=1 i=l 
N-l 

E ViSi = 0 


2 = 1 

z/fe = 0. 


( 180 ) 


48 










Since k was arbitrary, we have shown that Vj = 0 Vj. Therefore, since Ui is arbitrary, we 
have proven that is a basis. Finally, to agree with fll 75 p we choose 


z/. 




Cu'^Sj = 0 


( 181 ) 


(no sum on j), and it is clear that we have diagonalised this system to 


~ fxs 


Vj, 


( 182 ) 


at least up to a factor of an unknown positive coefficient C, which may in general depend 
on j. Thus, at least for the purposes of this proof, the system is decoupled. 

For ease of notation we dehne tq = for black holes and tq = 0 for solitons. We recall 
that ^j(ro) = 0 and ^'(ro) 7^ 0 in general. Therefore, noting that C, r^, /r, S and uf are 
all positive, then by integrating both sides it is easy to establish that for r > tq, S^'(r) 
has the same sign as i’j(r) - i.e. S^j(r) is monotonic; increasing (decreasing) if ^j(ro) > 0 

(q(ro) < 0).n 


Comment We briefly note that we unfortunately cannot transfer this result into telling us 
something similar about aj(r), since the are the differences between each pair of successive 
functions aj{r). In the case of su(2) alone, there is only one Sj = S and only one independent 
EGF oi = a, and since Oi + 02 = 0 for su(2), we hnd £ = ai — a 2 = 2ai = —2a2 = 2a. 
What is more, a is monotonic in the su(2) case^. Therefore, both aj{r) are monotonic 
there, and hence for embedded solutions, Sj = 2a are monotonic Vj. So this result does 
imply at least that any su(2) embedded solutions will also have all aj monotonic; though we 
note that due to tracelessness and the transform itself, various aj will be positive (negative) 
and monotonically increasing (decreasing) for r > for black holes or r > 0 for solitons, 
with at least one aj of each sign. 


A. Global regularity for p(r) > 0 

We now need to prove that the local solutions that we found will remain regular if we 
begin with their initial conditions and integrate out arbitrarily far. We again let tq = Vh 
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for black holes, and ro = 0 for solitons, so that we are discussing regularity in the range 
R = (ro,cxo). Also, note that we require the metric function /i(r) > 0 for all r E R; but 
we expect R to correspond to the external region of a black hole (and /i(r) > 0 Vr G i? for 
solitons), therefore taking p(r) > 0 Vr G i? is a natural physical requirement. We state the 
proposition: 


Proposition 16 As long as nir) > 0, the field equations / fldl \TU. \2l\ f^) remain regular 
in all field variables throughout the range r^ < r < oo for black holes, and 0 < r < cxd for 
solitons. 


Proof Take some ri > ro, and define intervals Q = (ro,ri) and Q = (ro,ri]. Since we 
have proven existence in some neighbourhood of the boundary points represented by r = ro, 
our strategy is to assume that the field variables remain regular on Q, and use the held 
equations to prove that they continue to remain regular on Q (i.e. at r = ri), as long as 
/r(r) > 0 Vr G Q, and thus we may integrate the held equations out arbitrarily far and into 
the asymptotic regime. 

We hrst note that all the terms on the right-hand sides of the Einstein equations flT^ and 
(fT6|) are non-negative, and therefore all of those terms are bounded below by 0 and above 
by the left-hand sides. This also implies that for non-trivial solutions (with non-constant 
m{r) and *S'(r)), m'{r) > 0 i.e. that m(r) is monotonically increasing (since m{ro) > 0); and 
that S'S~^ > 0 i.e. S"(r)S'(r) > 0, and since we can hx Sq to be positive, then S'(r) is also 
monotonically increasing. Finally, an important upper bound is given by fi{ri) > 0, which 
implies 

3 

2m(ri) < kri + (183) 

The right-hand side of this is positive (recalling the minimum event horizon radius for k = —1 
61), SO that m(r) is regular on Q and therefore so is /i(r). Now we consider (IT^ : due to 
the preceding comments, we may write 

2m'(r)>2/iG + ^. (184) 
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Now /i(r) must have a minimum value on Q, so define 


/imm = min{/i(r) : r e Qj. 


( 186 ) 


Then, integrating both sides of (11841) . we obtain 


ri 


2[m(r0 - m(r„)| ^ 


Mr. 


2/^252 


(186) 


ro 


and direct integration of (fT6|) gives us that In |S'(r)| and therefore S(r) is bounded in Q. 


Therefore, each of the integrals on the right-hand side of (11861) must also be bounded. 
So, using the Cauchy-Schwartz inequality: 



(187) 


and since the left-hand side is a sum of positive terms and the right-hand side is bounded 
above by fll86p . we can say that ujj(r) is bounded on Q for all j G {1,..., iV — 1}. 


In a similar manner we directly integrate (IT^ and consider only the first term on the 
right-hand side, which in turn is bounded above by the bounds on m{r): 


[m(ri) — m(ro)] > 


4^2 


af{r)dr. 


(188) 


ro 


j=l 


Extracting factors of bounded (positive) functions, it is elementary to show that 


v-i 

,.'2/ 


[m(ri) - m(ro)] > j y ] af{r)dr. 

i=i 


(189) 


0 


ro 
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However, again using Cauchy-Schwartz, we find 


N 


N 


i=i i=i 



( 190 ) 


and so likewise aj{r) is regular on Q for all j G {1, which in turn means that so is 

Sj(r), by definition. 


Finally we examine the Yang-Mills equations in the forms 022 p and 0291) . We begin by 
rewriting them as: 

/ ^ /\/ —SW^oji I ^ 

(192) 

Begin with fll92p . Note that all terms on the right-hand side are bounded on Q, so we can 
immediately write 

ri 

lx{n)S{n)u'j{n) = /i(ro)^(ro)a;'(ro) - j (193) 

ro 

and therefore the left-hand side is bounded, and so ca'(r) is regular on Q for all j G 

{1,...,Y-1}. 


Finally, taking fll9ip . the same is true: we may write 


2 2 
-Mri) = + 


ri 


S{n) 


S{r,) 


1 

~i^S 


Zjdr, 


(194) 


ro 

showing that since the right-hand side is bounded, the left-hand side will be too. Therefore 
£’j(r) is also bounded on Q for all j G {1, ...,N — 1}. Thus Proposition [Tbl is proven.□ 
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B. Asymptotic behaviour 


However, our preparations are still not quite complete. A big difference between the A = 0 
and A < 0 EYM cases is the dearth of solutions in the former case and the abundance in the 
latter, and this is entirely due to the different geometry of the manifold in the limit r —)■ cxd. 
Therefore, it is important to examine the behaviour of the equations in the asymptotic limit. 
We begin with the Yang-Mills equations in the same form as in Section IIV Al 





-SW.U. 1 
-^ 

4/iA ^ 


(195) 

(196) 


2 

First we let r —)■ cx), so that /r(r) —)■ ^ and S{r) —>■ 1, and then we make the co-ordinate 
change r = to convert the asymptotic field equations into autonomous form. This 

yields 


CpUj 

dr'^ 


-C (w--, + u,. 


(197) 


It may be seen that although these equations are clearly autonomous, they are not scale- 
invariant as they contain a reference to £; however this is not a problem as we are not 
interested in the limit £ —)■ cx) (i.e. |A| —)■ 0), and the equations remain well-behaved if we 
take ^ arbitrarily small (which we shall consider later). Suffice to say that equations (11971) 
are easily solved to hnd the critical points ujj and £j - the relevant equations are 


(198) 

(199) 

and certain elementary solutions are given as follows: 

(i) For all values of fc, we get the solution ujj = 0, £j arbitrary; 

(ii) For k = 1 only, we get the solution Qj = j{N — j), £j = 0; 
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where (i) is a metastable centre in the phase space and (ii) are saddle points, as we may 


expect^. Note that these solutions are independent of the value of i. Because (119811199p 


are non-linear equations, there may potentially be other critical points; but it can at least 
be observed that for each critical point, the given value for ojj solves the system (119811199p 
if and only if the given value for £j solves (119811199p . 

We notice that the study of the solutions locally as r —)■ oo implied no such constraints on 
the parameters. However, this is due to our choice of autonomous parameter, where r oc r~^. 
Similar to the purely magnetic adS cases^i^, this means that if we consider our solution as 
a trajectory in the phase space of the 4A^ — 4-dimensional system 



( 200 ) 


and then consider the transform to the variable r, then we can see that integration over 
the range r G [^ 1 , 00 ) becomes integration over r G [0,ri], and if ri is large enough to be 
considered asymptotic then the corresponding trajectory in terms of r will be very short. 
Hence the solution will not move all the way along the trajectory, and therefore the values 
of the gauge held at inhnity are unconstrained. We can compare this to the hat space case, 
in which there is a very diherent situation: the parameter used there must be such that 
r oc Inr, hence (ri,oo) 1 —)■ (ri,oo), and so in that case every solution must proceed to the 
end of its trajectory. This is the reason for the scarcity of solutions in the case A = 0, where 
in the purely magnetic case it is proven that solutions can only be found for certain discrete 
values of the gauge held parameters at inhnity. 

V. GLOBAL EXISTENCE OF NON-TRIVIAL SOLUTIONS 

Now we have all the machinery we need to construct global solutions which are genuinely 
non-trivial, i.e. solutions which do not appear in our list in Section fllBI In this penultimate 
section, we present arguments that allow us to piece together the local solutions we found 
at the boundaries in Section m using the global regularity results in Sections HVl and thus 
describe global non-trivial solutions to the held equations (1T5| fT6| [22il . We consider 
existence in several regimes, including in the limit of |A| —)■ 00 . 
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A. Existence of non-trivial solutions near existing solutions 


The crux of the proof of global existence is in the following powerful proposition, the 
essence of which is the proof that dyonic solutions to the field equations exist in open sets 
of the parameter space: 

Proposition 17 Assume we have an existing solution to the dyonic su(iV) topological field 
eguations where each MGF ujj{r) has Rj nodes and gauge fields have initial 

values {oojh, £'jh} for topological black holes or {/3fc, 9^} for spherical solitons. Then all initial 
gauge field values in some sufficiently small neighbourhood of the existing solutions will also 
give an su(A) dyonic topological solution to the field equations (a spherical solution for 
solitons) in which each MGF oj fir) has Rj nodes. 


Proof Assume we know of an existing solution to the dyonic 5u{N) topological field equa¬ 
tions, where each Ujfi) has Rj nodes. For black holes we have initial conditions {ujh fi 
0,Sfi fi 0}, and for solitons, {ffi fi 0,9k fi 0} (in the general case). From these initial con¬ 
ditions, Proposition [16] and Section IIVBI it show that as long as /i(r) > 0 we may integrate 
this solution out arbitrarily far into the asymptotic regime to obtain a solution which will 
satisfy the boundary conditions as r —)■ cxd. For the rest of the argument, we assume that 
i and rh are fixed (where Vh = 0 for solitons); that each MGF Uj has Rj nodes; and that 
again, ro = rh for black holes and ro = 0 for solitons. 

From the local existence results (Propositions [HI [IS! and [TTD . we know that for any set 
of initial values there are solutions locally near the event horizon (for a black hole, or the 
origin for a soliton), and that solutions are analytic in their choice of initial conditions. For 
an existing dyonic su(A) solution, it must be true that /i(r) > 0 for all r G [ro, cxd). So, by 
analyticity, the nearby dyonic 5u{N) solution will also have /i(r) > 0 for all r G [ro,rc] for 
some r = rc with vq Kr^ < oo. By Proposition [T6| this nearby solution will also be regular 
on [ro,rc]. 

Now, choose some ri >> ro such that for the existing solution, m(ri)/ri << 1. Let 
{(jjjh,£'jh} ^ black hole solution, or {fik,9k} for a soliton) be a different set of initial 
conditions at r = ro in some neighbourhood of the existing solution, and let m(r) be the 
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mass function of that solution. By analyticity (as above), these will also be regular on [ro, ri] 
- i.e. /i(r) > 0 on this interval - and the MGFs Uj will each have Rj nodes. 

Also it is then the case that rh{ri)/ri << 1, and since ri >> tq we consider this the 
asymptotic regime. Provided vi is large enough (and hence ti is very small), the solution 
will not move very far along its phase plane trajectory as ri —)■ oo. Therefore m{r)/r remains 
small and the asymptotic regime remains valid. Therefore the solution will remain regular. □ 

Comment We note that the set of generally dyonic solutions trivially includes the set of 
purely magnetic solutions, simply by choosing those ones for which £j = aj = 0 Vj. Given 
then that this result also applies to existing topological solutions, su(2) dyonic solutions and 
purely magnetic su(iV) solutions, this is a really rather general result that we can use to 
prove existence in a number of regimes, using these established solutions^i^i^^i^. We also 
recall that we are most interested in solutions for which the MGFs have no nodes, so this is 
the kind of solution we look for here. 

We sum up this subsection in the following theorem, which uses Proposition [T7] and known 
trivial solutions and solutions from previous research to generate nearby non-trivial solutions. 
A similar structure of argument will apply in each case, so we choose to use a shorthand 
for brevity and clarity: X ^ Y will be shorthand for “the existence of solution X implies 
the existence of the nearby solution Y using Proposition [T7I’ • (For a longhand version of 
the argument structure, see Theorem 8 in our previous work—). Note that our Proposition 
[m relies on solutions being analytic in some neighbourhood of existing solutions, so while 
that guarantees that nodeless solutions will have nearby nodeless solutions, it also means 
that if we use a purely magnetic solution as our ‘trivial’ solution, we will get a neighbouring 
solution for which all the £j (and hence due to their zero sum) will be smalt, and if we 
use an su(2) embedded solution we get all aj monotonic. We have tried to highlight these 
properties in each case for the sake of completeness. 

Theorem 18 Regimes of existence of non-trivial nodeless solutions 

Non-trivial global solutions exist to the field eguations 1751 ITR. \2l\ 17^) in the following 
list of regimes found nearby existing solutions. 

• For strictly topological black holes (i.e. with k 1): 
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1. Topological 5u(2) purely magnetic solutions^, all of which are nodeless 

Nearby non-trivial su(2) topological dyonic solutions, all of which are nodeless 
due to analyticity, and for which the only independent EGF a is monotonic and 
small - and hence embedded su(iV) topological dyonic solutions based on these 
(by dUD) 

Non-trivial (i.e. non-embedded) solutions nearby embedded 5u(iV) topological 
dyonic solutions, again all of which are nodeless due to analyticity, and for which 
all ttj are monotonic and small. 

2. Topological su(A^) purely magnetic solutions^, some of which are nodeless 

Nearby non-trivial dyonic 5u(iV) solutions, some of which are nodeless and for 
which all aj are small, and not necessarily monotonic. 

• For spherically symmetric black holes and solitons (i.e. for k = 1): 

3. Non-trivial dyonic su(2) solutions^, some of which are nodeless, imply the exis¬ 
tence of 5 u(A^) embedded dyonic solutions (HTil 

Non-trivial (i.e. non-embedded) solutions nearby su(iV) embedded dyonic 
solutions, some of which are nodeless, and for which aj are all monotonic but 
otherwise quite general. 

4. Non-trivial purely magnetic su(A^) solutions^S^ some of which are nodeless 

Nearby non-trivial su(A^) dyonic solutions, some of which are nodeless, and 
for which all aj are small, and not necessarily monotonic. 

5. For black holes only: The su(A^) SadS purely magnetic solution (see Section ing , 
which is nodeless 

Nearby non-trivial solutions all of which are nodeless, and for which all aj are 
small and not necessarily monotonic. 

6. For solitons only: The su(iV) pure adS solution (see Section Hi BD . which is node¬ 
less 

Nearby non-trivial su(iV) dyonic solutions, all of which are nodeless, and for 
which all aj are small and not necessarily monotonic. 
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B. Existence of solutions as |A| ^ oo 

So far, we have proven the global existence of non-trivial solutions for any hxed value 
of the cosmological constant A < 0. However, we see that numerical results in the purely 
magnetic case^*^ imply that for hxed A, the gauge held parameter space in which we hnd 
nodeless MGFs shrinks as N increases. However, we also hnd that if we make |A| large 
enough, all solutions we hnd have nodeless MGFs. In addition to this, we note that in the 
case of purely magnetic su(2)^ and 5u{N)^ solutions, |A| —>■ oo emerges analytically as a 
condition of stability in the gravitational sector. 

We have not been able to prove that all solutions for |A| —)■ cxo are nodeless, but we could 
at least prove nodeless solutions exist in this limit^. Since the purely magnetic solutions 
are a limiting case of our dyonic solutions, we expect these results to carry over to at least 
some dyonic solutions, i.e. those with the EGFs all small. Hence we do not attempt to 
prove that all solutions are nodeless in the MGFs as |A| ^ oo] rather, we attempt to prove 
that we can hnd at least some such solutions, i.e. that this a sufficient condition for the 
existence of nodeless solutions. We hnally note that the shrinking parameter space for hxed 
A and increasing N implies a possible necessary condition for nodeless existence involving 
the initial gauge held parameters, though this is something we have not yet been able to 
identify. 

Therefore, motivated by the above results and by analogous analytical results for existence 
in this regime^i^, we now describe transforms of our held equations that will allow us to 
take the limit |A| —)■ oo (i.e. i ^ 0) sensibly. We emphasise that we cannot let £ = 0 as then 
the asymptotic parameter r becomes non-regular. Here, it is more convenient to take black 
hole and soliton solutions separately. 

Proposition 19 There exist non-trivial dyonic solutions to the field equations mmm 
analytic in some neighbourhood of i = 0, for any choice of initial MGF values ojjh (or 
fik) and for initial EGF values in some neighbourhood of = 0 Vj (or 6k = 0 Vfcj. 


Proof We begin by considering black holes. In the purely magnetic topological case, we 
were able to deduce the existence and uniqueness of black hole solutions as £ —)■ 0, by using 


the following rescaling to simplify the equations: 


m = mi‘^, fi = iiP = + r^. (201) 

r 

Unfortunately, in our case this does not simplify the equations much at all; all we can prove 
is that m must once again be a constant. However, we know that for the su(iV) purely 
magnetic equations, we obtain the following unique solution (for all r > r/j): 

rri(r) = Ir)), S{r) = 1, Uj{r) = Ujh- (202) 

This is valid for all initial MGF values Ujh- Therefore, if we append Sj{r) = 0 to this 
list, we will end up with a (trivially) dyonic solution as £ ^ 0. (Due to non-linearity it is 
possible that there are other solutions for Sj{r) ^ 0, i.e. that this solution is not unique, 
however we at least see that fl202p solves the transformed equations if and only ii £j = Q.) 
We recall that £j = 0 Vr plus the zero-sum of the aj implies that aj = 0 Vr also. 

It is then also clear that the argument in Proposition 11 in previous work— will carry 
identically across for this solution, thus we have a ‘trivial’ solution to work with for any 

including ^ arbitrarily small. Hence we shall take this solution and show that by hxing 
rh and oojh, letting i be arbitrarily small and varying only we can hnd solutions for 
non-identically-zero EGFs. 

It is clear from the dehnition of /i that i only makes a difference as r becomes large, 
and is ignorable otherwise. Therefore, if we can prove local existence for £ ^ 0 as r ^ oo 
by adapting Proposition [HI then Proposition [T7] can be used to prove nearby non-trivial 
solutions exist in some neighbourhood of these purely magnetic solutions (i.e. for £j, i small 
but non-zero). 

We apply (120Ih to Proposition [HI and dehne the following new quantities: 

Ai = eXi, Woo = Woo = (203) 

Examining the proof of Proposition [HI we see that all of the polynomials ^ 00 , 2 , "Hoopj, ^ 00 , 2 j 
change internally but remain regular (including as £ —)■ 0). We must also replace the equation 
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x{dA/dx) = 0 with 



(204) 


but the only equation which changes is that of Ai, which becomes 



(205) 


Therefore we see that this proposition remains valid as £ —0. As for the asymptotic regime, 
we can see that our argument in Section IIVBI remains valid for all i small but non-zero. 

The argument that proves the existence of black hole solutions in the regime is essentially 
identical to that in Proposition [171 We £x Vh and uj^h non-zero and we £x i arbitrarily small. 
We then choose some ri >> so that we can consider ri in the asymptotic regime, and 
we consider varying only Sj Then Propositions IT^ and fTbl conhrm that for Sj sufficiently 
small we can hnd solutions near purely magnetic solutions which will begin near r = and 
remain regular in the range (r/j,ri], and that those solutions will be nodeless in the MGFs 
due to analyticity. Finally, once we are in the asymptotic regime, we can use the logic in 
Section IIVBI to ensure that solutions will remain regular as r —>■ cxo and that all ujj will have 


no nodes. 


Luckily in the case of solitons, the argument is very similar. The analogous argument 
for existence of purely magnetic solutions in this regime exists already^^, so we will just 
mention that we must be more careful about how we take the limit £ —)■ 0 since we must 


use the parameter r~^ at inhnity. For black holes this is hne since > 0 and therefore r~^ 
is bounded and thus we have a natural ‘scale’ to work with, but for solitons, we must have 
r = 0 at some point so that r~^ is not bounded. 


There, we dehned the following transforms. 


r = ix, m{r) = ifh{x), 


(206) 


and found that for the MGFs, it was easier to work with functions dehned by the basis of 
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eigenvectors that we used near the origin: 


N N 

cu = LJo + X] Wk{r)vkr^ = cjq + ^ Wk{B)tx^Vk. (207) 

k=2 k=2 

In that case, the solution turns out to be 

m(r) = 0, ^(r) = 1, Wk (x 2 F 1 —-—;. (208) 

The functions vuk are essentially the same as (3k but for a factor of jj. It can be noted that 
the properties of hypergeometric functions can be used to establish that (12081) also fulhls the 
boundary conditions at r —)• cxo. Using a very similar argument to the black hole case, we 
append £j = 0 Vj, i.e. 9k{x) = 0 V/c, to our solution (I208p : and we deduce that we can hnd 
global solutions regular and analytic in all held variables for i arbitrarily small, for arbitrary 
values of ^k and for 9^ small (VA;).n 


VI. CONCLUSIONS 

The purpose of this research was to investigate the existence of black hole and soliton 
solutions to topological adS dyonic equations in four-dimensional su(iV) EYM theory, moti¬ 
vated by previous existence results for purely magnetic spherical su{N) solutions'^S, purely 
magnetic topological su(V) black hole solutions^E^ and dyonic su(2) spherical solutions^. 

We began by using a previously derived gauge potential appropriate to the case^. We used 
this to derive the held equations in this case - two Einstein equations and 2N — 2 independent 
Yang-Mills equations - and found several trivial solutions, including an embedding of su(2) 
in su{N) (Proposition [1]). At r = and r —)■ 00 , we used physicality requirements to 
establish appropriate boundary conditions in an elementary fashion. For solutions regular 
at the origin, the situation was much more complicated and required detailed analysis. 

Upon expanding the equations in power series at the origin, we discovered that to establish 
consistency of solutions, we needed to solve a tri-diagonal system similar to one previously 
considered by Kiinzle^. In our case we have a mutual system of equations, but we found 
that each set may be expanded in the same basis of eigenvectors, which simplihed things a 
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little. We used a form of ‘ladder’ induction over two induction sentences to prove that each 
expansion parameter in each system depended only on previous parameters, and that for 
each gauge field we needed — 1 independent parameters to describe the power series at 
r = 0 (Proposition E]). 

We then proceeded to prove the existence of solutions locally near the boundaries (Propo¬ 
sitions m Us] and [14]), which are regular and analytic in their boundary values, using an 
established theorem of ordinary differential equations (Theorem [2|). Existence is elementary 
to establish in the cases oi r = and r —)■ oo, and a little more complicated for regular 
solutions at r = 0. Nonetheless we establish local existence here too. 

After that, we constructed a series of arguments which form the heart of the proof. After 
proving that all £j are monotonic (Proposition IT5l) we proved that if a solution remains 
regular over a small interval near the event horizon (origin), then with the condition /i(r) > 0 
for all r > r/j (r > 0) we may continue to integrate that solution out regularly into the 
asymptotic regime fProposition [T6|) . We examined this regime and discovered that due to 
the parameter we used to render the asymptotic equations autonomous, the solution would 
continue to remain regular for r arbitrarily large (Section IIVBI) . Finally, we proved that 
for these held equations, solutions exist in open sets (Proposition [T7]); and given that the 
literature is now fairly abundant with known solutions that result from trivialising these 
held equations in some way^'2^‘^>^, we are able to prove the global existence of non-trivial 
solutions in a number of regimes (Theorem IT8|l . In addition, we proved that solutions 
could be found in the limit |A| —>■ cxo for arbitrary initial parameters for ujj, and for initial 
parameters of £j small (Proposition . 

The main result of this paper is the proof of existence of nodeless non-trivial black hole 
and soliton solutions to four-dimensional topological su(A^) dyonic EYM equations in various 
regimes: nearby existing solutions, and in the limit |A| large. In particular, we have shown 
that we may dress a black hole or soliton with an arbitrarily large amount of gauge held 
hair. There are several future research directions that are suggested by the results here. The 
possibility of dressing a dyonic black hole with arbitrary amounts of hair suggests there may 
be some work to do in extending the gauge group as large as we possibly can, and therefore 
the dyonic su(oo) case becomes interesting. We note that the purely magnetic case has 
already been considered^, in which evidence of the existence of solutions is provided. String 
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theories are characterised by enormous gauge symmetries and hence provide a motivation 
for this work. 

Another possible extension we could make is suggested by Gubser— , who found some 
very interesting results for dyonic su(2) planar black holes in adS; namely, a second-order 
phase transition between the embedded planar RNTadS black hole and a black hole with a 
non-trivial Yang-Mills held condensate. It would be natural to ask how his results generalise 
to su{N) in light of the solutions that we have discovered here. 

In addition to this there are questions which arise about the impact of black hole hair 
on other areas of gravitational physics. For instance, it may be valuable to consider the 
adS/CFT (Conformal Field Theory) correspondence in light of this work—'—, since it has 
been conjectured that there are observables in the dual CFT which are sensitive to the 
presence of black hole hair— It would also be of interest to know whether these topological 
models would be valuable to modelling holographic superconductors, since planar black hole 
models have recently been used in this research areai^. 

Finally, there is the important question of further conhrming or rehning the generalised 
‘No-hair’ theorem, on which this work will have direct implications. The statement of this 
theorem given by Bizort^ is: 

“In any given matter model, stable black holes will be characterised by a hnite 
number of global charges.” 

Therefore, the next important problem here is to examine the stability of the solutions that 
we have found, since we can establish the existence of stable solutions in the purely magnetic 
spherical case^, and results have just emerged proving the stability of some topological 
black hole solutions^ and su(2) dyonic solutions^. We expect however that establishing 
the stability of our dyonic solutions will be a highly non-trivial problem, since the presence 
of the electric gauge held will disallow the decoupling that happens in the purely magnetic 
case, making it much harder to prove analytical assertions. 
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